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Exact spectra of periodic samples are computed up to A^ = 36. Evidence of an 

> ■ 

QQ ■ extensive set of low lying levels, lower than the softest magnons, is exhibited. 
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f^ , These low lying quantum states are degenerated in the thermodynamic limit; 

^D I their symmetries and dynamics as well as their finite-size scaling are strong 

^I* ■ arguments in favor of Neel order. It is shown that the Neel order parameter 
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agrees with first-order spin- wave calculations. A simple explanation of the 
low energy dynamics is given as well as the numerical determinations of the 
energies, order parameter and spin susceptibilities of the studied samples. It 
is shown how suitable boundary conditions, which do not frustrate Neel order, 



k> I allow the study of samples with A^ = 3p -|- 1 spins. A thorough study of these 

C^ I situations is done in parallel with the more conventional case N = 3p. 
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I. INTRODUCTION 

The nature of the thermodynamic ground-state of the spin-1/2 Heisenberg antiferro- 
magnetic hamiltonian, in two space dimensions, is still an open question. There have been 
considerable amounts of theoretical or numerical works on Heisenberg antiferromagnets[l-45] 
but few exacts results. It is known that in one and two dimensions the system is disordered 
at T 7^ III], and that the one-dimensional system does not exhibit Neel order even at 
T = 0. The two-dimensional case is more controversial. There is a rather large consensus 
on the existence of Neel order at T = on the (unfrustrated) square lattice P-p!2|. The 



situation is much more puzzling as regards the triangular lattice case. It was indeed the 



first system to be proposed by Anderson [T^ and Fazekas |]I4[ as a candidate for a spin- 
liquid. On this lattice the "frustration" implies that the classical system is not very stable 
{Ed =< 2si.Sj >= —1/4), and the spin-wave calculations predict an important reduction 



(by about one half) of the sublattice magnetization by quantum fluctuations |T5|-|T^]. Per- 



turbation theory [ffH], series expansions [19| and high temperature calculations pO[ have 



been developed which suggest that the spin-wave calculations possibly underestimate this 
renormalization. Many variational calculations have been done exhibiting either ordered 



,2T| or disordered solutions p2|-p4 



In the square lattice case, numerical methods (Q.M.C, Ulam's or Trotter-Suzuki meth- 
ods) have brought very interesting indications on Neel order 0,^11]. Unfortunately these 
methods which allow to handle large samples cannot be applied to the triangular case: they 
lay on a property of positivity of offdiagonal matrix elements which is violated in the triangu- 
lar case; it is the well known sign problem which plagues many studies of strongly correlated 
fermions. Exact diagonalizations of the hamiltonian are thus the last resort to gather new 



information on these models. This approach has been developed by other authors p5H32 
Most of them conclude to the absence of Neel order for the triangular Heisenberg antiferro- 
magnet (THA) at T = 0. But two requirements are to be met to analyze the raw numerical 
data: a consistent finite-size scaling analysis and a quantum definition of observables. With 
these two constraints we will show in this paper that all the numerical data point to an 
ordered ground-state for the THA. 

As the second and perhaps more important objective we want to illustrate the spectral 



properties of an ordered quantum antiferromagnet on finite lattices and how they embody 
the characteristics of the symmetry breaking state (parts of these results have already been 
pubhshed H). 

In section II, we first show the numerical spectra of the THA on periodic samples with 
N=9 up to N=36 spins and we present the essential characteristics of these spectra. Their 
low lying levels form two families: the first one contains an extensive number of states 
called QDJS in ref |3^, which collapse to the ground-state as A^~^, and are characterized 



by the spatial symmetries of the classical Neel ground-state. The second family collapse 
more slowly to the ground-state as N~^^'^ and forms the familiar one magnon excitations. 
Such a structure has been conjectured a long time ago by Anderson in his seminal paper 
on antiferromagnets 0, and the subject has already been studied in the square lattice case 
J3ff7|. 



Section III is devoted to the understanding of the QDJS of the triangular Heisenberg 
hamiltonian. We show how the hypothesis of a Neel order explains the number of levels 
appearing for each spin, their symmetries and their dynamics. The analysis of these spectra 
gives a first information on the spin susceptibilities x± ^ind x\\ of the THA. In order to 
check our general ideas, we extend the analysis to samples with N = 3p + 1 spins. Suitable 
boundary conditions allow us to study situations which do not frustrate Neel order. We 
compare all the results of energies and susceptibilities with the first-order spin-wave results. 

In section IV, we analyze the static spin-spin structure function and the Neel order 
parameter. We show how diagonalization results sustain the hypothesis of Neel order and 
how the analysis of some previous authors was in error, leading to confusing conclusions. 

We conclude by a brief discussion of the interest and limitations of this approach. 

In the Appendices, we develop the group theoretical analysis of the THA and explain 
our numerical method: this is a technical but crucial point in order to obtain the needed 
information in a minimum computer time, within memory capacities of today computers. 
We had for example to make sure that the low lying levels of all the different irreducible 
representations order as they should in the framework of our hypothesis. This leads us to 
compute the whole spectrum of the samples up to A^ = 21 and a large number of low lying 
levels in each irreducible representations for all sample sizes up to A^ = 27. 



II. SPECTRA OF PERIODIC SAMPLES: AN OVERVIEW 

All the results presented here have been obtained by diagonalization of the spin-1/2 
Heisenberg hamiltonian of periodic samples on the triangular lattice ( p5[ and Appendix A 
for more details). The Heisenberg hamiltonian reads: 

^ = E 2 ^^■^J (1) 

where the sum runs over pairs of nearest neighbors and Sj, Sj are the spin-1/2 operators on 
sites i and j. 

In this section, we are concerned with periodic samples with N = 3p. Such samples do not 
frustrate the classical Neel order. A classical Neel state on the triangular lattice has coplanar 
spins with a three-fold rotational symmetry defining three magnetic sublattices {A, B, C) 
(see fig.|T]). On each sublattice, the spins are ferromagnetically aligned and the angles between 
the magnetizations of two sublattices are ±2tt/3. The total spin of a triangular plaquette is 
zero and the rotations of the spins around an ABC triangle could be clockwise or counter- 
clockwise defining two opposite hehcities. The point symmetry group of these classical 
solutions is Cs^, and the translational symmetry group is that of a sublattice, say A. These 
classical solutions break the translation and spin rotation invariance of Eq.|l] whereas the 
quantum eigenstates that we will now consider do not. 

Using group theory, described in Appendix A, we have computed the complete spectrum 
of Eq.|l] for A^ = 9, 12, 21. The 5 lowest energies have been computed for A^ = 27 in all the 
irreducible representations (IR), and the 3 lowest energies for A^ = 36 in the homogeneous 
states k = 0, invariant under rotation. The energy spectra are given in fig.|^ and the lowest 
levels listed in table |. Our results are in perfect agreement with previous diagonalization 
results P^g. 



A first very striking feature of these spectra could be read in fig.|^: the lowest energy 
levels in each spectra order with increasing S. This is strongly reminiscent of the theorem by 



Lieb and Mattis for bipartite lattices PS|. But there is presently no proof that the property 
holds for other lattices. We find this is true for other non bipartite problems like the Kagome 
Heisenberg antiferromagnet or for the Ji — J2 models on triangular and Kagome lattices. 



The second striking feature of these spectra is the existence in each 5" subspace of a 
family of low lying levels well separated from the others: we called these states the "Quasi- 
Degenerate- Joint States" (QDJS) in ref. |Q. In fig.||, it is shown that the QDJS energies 
(6A^ < Sj.Sj >) stand around a line E = S{S + l)/(2/iv), where the moment of inertia In 
is essentially proportional to A^ (see fig.^)- This suggests that these levels collapse as A^~^ 
to the ground-state. We shall show in the next section that this family of 0{N^) QDJS has 
all the properties expected for the description of Neel quantum ground-states: in particular, 
the Csv and magnetic sublattice translation invariance of the Neel states (they only contain 
the k = or ±ko IRs of the lattice translation operators where ±ko are the two wavevectors 
of the corners of the crystal Brillouin zone, mapping on the center of the magnetic Brillouin 
zone). The symmetry breaking Neel states are linear combinations of these QDJS. 

To study higher excited states, let us first come back to the A^ = 9 spectrum of fig.§. 
Above the basal line of QDJS one sees very clearly two families of levels. 

The first excited-state family (horizontal bars in fig.^) consists of eigenstates of the 
translation operator with k 7^ and k 7^ ±ko, they are typically states involved in one 
magnon excitations. What is usually called a magnon in solid state physics is in fact a spatial 
modulation of the Neel state. As the Neel states should be seen as linear combinations of 
the QDJS, the one magnon excitations are in fact linear superpositions of this first family 
of excited states. A one magnon excitation is thus pictured in the spectrum of fig.^ as a 
AS = 1 collective excitation of the QDJS (three legs symbols family) towards the k family 
of levels (this picture has been ascertained by the computation of the dynamic structure 
factor of the first QDJS of the A^ = 12 sample |3|). 

The second family of excited states (black triangles) belongs to the non trivial IRs of C3 
[7?.27r\E' = exp(±z27r/3)\E'], this is the first family of levels (labeled 3— in the following) that 
could sustain chiral states (7^27r/3 is a spatial rotation of angle 27r/3). The last family of 
excited states (open triangles) belongs to the non trivial IR of C^y, the states invariant under 
a rotation of 27r/3 and odd under axial symmetry {a^) [R.221'^ = \I', cra;\E' = — \I'] (labeled 
3 = X—, see table III). 

For larger A^, the spectra become more and more dense, but this hierarchy of levels is 
always obeyed. The magnon states collapse to the ground-state roughly as k i.e. as N~^^'^ 



whereas the energy of the chiral states seems to have a gap {E^^ — Eq = 6, 6, 1.908, 2.958, 
0.7661, 2.454, 2.641, 2.01, 2.546, for A^ = 7, 9, 12, 13, 16, 19, 21, 25, 27). The hypothesis 
of a chiral ground-state on the triangular lattice seems excluded by our results [^pl . The 
energies of the 3 = x— states are much higher in the spectra (out of the range of energies 
of table I). 

III. QUALITATIVE ANALYSIS OF THE QUASI-DEGENERATE 
GROUND-STATE MULTIPLICITY 

The symmetries and the dynamics of the QDJS are essential to understand the nature 
of the order in the thermodynamic limit. Let us begin our investigation by the symmetry 
analysis of the eigenspectrum. 

A. The QDJS as the coupling of three spins 

It is straightforward to verify (see for example Table |) that the QDJS family exhibits 
the exact number of states expected from the coupling of three spins of length A^/6: 

A^5 = min(2^ + l,A^/2-5+l) (2) 

This number is readily obtained by noticing that the Hilbert space of two spins Sa = Sb = 
N/6 can be split in (A^/3 + 1) subspaces associated to the eigenvalues of their sum Sa+b 
(with 15^ — 5*^1 = < Sa+b ^ Sa + Sb = N/3). The coupling of Sa+b with Sc gives then 
a total spin S {\Sa+b — N/6\ < S < Sa+b + N/6), which leads to result Eq.||. 

Such a property is not unexpected if we look at that family of levels as arising from 
the renormalization of the classical Neel state by quantum fluctuations. In this purpose 
it is most useful to Fourier analyze the hamiltonian Eq.|I]; separating the k = and ±ko 
contributions from the others, Eq.|l] reads: 

n = no + v (3) 

where 

7^o = |(S^-S^-S|-S^) (4) 



Sa (resp. Sb,Sc) is the total spin of the A (resp. B,C) sublattice, and 

V= Yl /(k)Sk.S-k (5) 

fc^O;±fco 



Sk = -y= ^ Si exp(ik.ri) (6) 



with 



/N 
and 

/(k) = exp(ik.Ui) + exp(2k.U2) + exp(— ik.(ui + U2)) + c.c (7) 

where Ui and U2 are the two basis vectors of the lattice. 

Ho commutes with Sa, S^, Sb, S^, S^, S^ and the eigenstates of Hq (written as 
I \l/o(?, 5, M5) >) are eigenstates of S^, S^, S^, S^- with eigenvalues: 

E{S, Sa, Sb, Sc) = ^ {S{S + 1) - Sa{Sa + 1) - Sb{Sb + 1) - Sc{Sc + 1)) • (8) 

For each 5* value the lowest eigenstates of Ho are obtained for Sa = Sb = Sc = ^^ their 
energies are: 

Eo{S) = ~^{N + Q) + ^S{S + l). (9) 

These states |\l'Q(i, 5*, M5) > (with i from 1 to A^^), fully polarized on each magnetic sub- 
lattice, are the projections of the classical Neel states on the various IRs of SU{2). They 
present the usual 25* + 1 degeneracy associated to Ms multiplied by the number Ns of 
different couplings of three spins A^/6 (Eq.|^). 

As S^, S^, S^ do not commute with V, |\I^o(i, S*, M^) > are not eigenstates of H, but 
we can look at them as the first (bad) approximation to the exact QDJS: the perturbation 
V dresses these "classical states" with quantum fluctuations, decreasing the average value 
of the sublattice magnetizations and lowering their energy towards the exact results. This 
process can lift the A^^ degeneracy of the |\I'o(i, S, Ms) > and it is indeed what is observed 
in the exact spectra of QDJS, but we verify in fig.^ and Table | that the set of QDJS has 
exactly the correct multiplicity Ns- 



B. Symmetries of the QDJS 

The QDJ eigenstates belong to the three following IRs of the space symmetry group: 

Ti : [ k = 0, 7^,^ = ^ , 7^2V3^ = ^ , a,^ = ^ ], T2 : [ k = 0, 7^,^ = -^ , 
7^27r/3^ = ^ , a^^ = ^ ] and Tg : [ k = ±ko , 7^2,/3^ = ^ , a^^ = ^ ], where 7^^ is a 
rotation of angle and o"^: is an axial symmetry. 

They appear with regular rules (described below in III.D and in Appendix B) in all the 
IRs of SU(2): that is for each S value. 

This proves that the QDJ eigenstates are invariant: 

1) under translations of the magnetic sublattices. They only contain the k = or ±ko 
IRs of the lattice translation operators. 

2) under the point group Cs^ (7^27r/3^ = ^, o-x^ = ^)- 

The appearance of the k = and ±ko IRs of the translation group in this quasi- 
degenerate ground-state multiplicity allows to build states which breaks the translational 
symmetry of the lattice (as it is the case of the Neel state). In the same approach, the 
appearance of all the IRs of C2 (states where TZn"^ = ±^) allows to break the inversion 
symmetry, whereas the presence of all IRs of SU{2) allows to build states with a vectorial 
magnetization pointing in a given direction. In summary, all the IRs which keep invariant 
the Neel states appear in the QDJS and no others. Note that the low lying levels of TYq 
described in the previous subsection belong to the same IRs than the QDJS. It is straight- 
forward to verify that the quantum perturbation V is also invariant under the symmetry 
group of the classical Neel states; thus the renormalization by V of the sublattices spins can 
take place without disturbance of the symmetries of the classical Neel states. 

The picture of the Neel states as linear combinations of the QDJS becomes plausible. 
But the symmetry argument does not prove by itself that the low lying levels of ?io have 
evolved identically in the dressing by the quantum fluctuations induced by V. For example, 
the simple question "do all the QDJS have the same extensive sublattice magnetization?" 
cannot be answered from only symmetry arguments. A still more basic issue is discussed in 
the following subsection. 



C. The effective dynamics of the QDJS 

We have shown in the previous subsections that the QDJS possess the same space sym- 
metries than the (sphericahzed) classical Neel states. We have now to make sure that the 
dynamics of this set of states can be reduced in the thermodynamic limit to that of a col- 
lective variable: the order parameter of the antiferromagnet. This is indeed a very serious 
issue and the necessary condition for the rigidity of the supposed-to-be ordered phase. 

To precise the nature of this order parameter, let us concentrate on a presumed Neel 
order on the triangular lattice. Defining a specific Neel state requires the knowledge of 
exactly three angles: two angles locate the helicity T defined as: 

T = Y^ (si A Sj + Sj A Sk + Sk A Si) (10) 

<i,j,k> 

where the sum is taken on upward triangles in the counter-clockwise direction, the third 
angle locates the direction of the magnetization of one sublattice. In this paper we shall call 
(3) the direction of the helicity, (1) the direction of the A-sublattice magnetization and (2) 
the third orthogonal direction. When the reference to a laboratory frame will be necessary, 
we shall add a prime to the Neel-axes frame indices, keeping the unprimed quantities for the 
laboratory frame. In the Neel ground-states (and in their first long-wavelength excitations) 
the length of the sublattices magnetization is supposed to be constant; thus, the orientations 
of the Neel frame are the only variables of the problem (homogeneous on the lattice in the 
ground-states and slowly spatially variable in the first excitations): in other words the order 
parameter is an element of 5*0(3). Let us now consider the dynamics of this collective 
variable in the homogeneous Neel state. On a finite lattice, this collective variable has a 
finite "inertia" and its free dynamics is entirely determined by its angular nature and the 
isotropy of spin space. We thus expect it to be that of a free top. Such a dynamics is at 
best described in the frame of principal axes of the object. The planar symmetry of the 
Neel state implies that one of the principal axes will be directed perpendicular to the plane 
of the spins, that is parallel to the helicity T. In the principal axes of the Neel state, the 
hamiltonian describing the free dynamics of the system reads: 

C2 c2 c2 

H^^ = ^ + ^ + ^ (11) 

'" 2/i 2/2 2/3 ^ ^ 



where 5*1, 5*2, S^ are the three components of the total spin of the system and Ii,l2, h the 
principal moments of inertia. Ii, I2, I3 are indeed hnear responses to homogeneous magnetic 
fields in the spin plane {Ii, I2) or perpendicular to it {I3). They are the extensive homo- 
geneous susceptibilities of this system. In the thermodynamic limit, we expect Ji and I2 
to be equal by symmetry, as they are associated to magnetic fields which rotate the spins 
out of the plane: we denote them I± = Nx±, and I3 is denoted I\\ = Nx\\, where x± and 
X\\ are the perpendicular and parallel susceptibilities. For the classical Heisenberg model, 
X± = X\\ = 1/18 (see Eq.^. Nevertheless, it is likely that on the triangular lattice the renor- 



malization of the two quantities by quantum fluctuations can be different fi^-^. Thus, we 
shall focus on the dynamics of a symmetric top. 

All the considerations up to now lay on the angular nature of the collective variable and 
are valid in a classical as well as in a quantum mechanics point of view. The quantization 
rules can be obtained from elementary quantum mechanics, but it is interesting to look first 
at the classical approach of the dynamics of this top [Q. The free dynamics of a top is 



entirely described by the relative motion of three directions (see fig.^: the total angular 
momentum S (a conserved quantity which is fixed in the lab frame), the angular rotation 
vector O and the principal axis of inertia of the top (axis 3 , directed along the helicity of the 
Neel state). For a free symmetric top, these three directions are in fact always coplanar, and 
the global motion is the combination of two rotations: the precession of axis 3 of the top 
around the total angular momentum S with an angular velocity Qpr = S/Ii and the uniform 
spinning of the top on itself around 3 with the angular velocity Q^' = S^' /I3 = S cosO/I^. 
The dynamics reduces to two separate motions of two angular variables constrained to vary 
on [0, Att] {SU(2) variables). In quantum mechanics, these two angular constraints imply the 
quantizations of the eigenstates of the associated rotation generators. The first condition 
provides us trivially with the quantization of the total angular momentum, and its (25* + 1) 
degeneracy. The second condition is not so trivial and implies the quantization of the 
projection of the total spin on the 3 axis of the top (not to be confused with the projection 
of the total spin on the 3 lab-frame axis). For a given 5* value, 5*3/ can thus take 25+1 
values ranging from S to —S. The hamiltonian then appears in the canonical form: 

10 



The degeneracy of the eigenlevels of Eq.|T2| is {2S+1) for S31 = and 2 x (25' + 1) for Sy 7^ 0; 
for a given S value of the total momentum the dimension of the Hilbert space is (25* + 1)^. 

Thus the rigidity of the ordered Neel states, implies that the low lying spectrum of 
Eq.|l| must map in the thermodynamic hmit, on the spectrum of Eq.|l2[ that is indeed a very 
striking feature of finite-size spectra displayed in fig.^ and |^ where both the leading behavior 
of Eq.|l2| and the global multiplicity (25* + 1)^ of the top is well verified up to a total spin 
S = N/6 (for higher values of S, the multiplicity (25* + l)Ns {Ns given by Eq.|^) is lower 
than {2S + 1)^). The inertia I± in Eq.|T^ is indeed an extensive quantity scaling as A^ (see 
fig.Q). This leads to a determination of Xi. = ^l/N which is not very different, and only a 
bit smaller than the first-order spin- wave calculations of Chubukov et al. |^ (see fig.|^). 

The situation is more difficult as regard the precise determination of the anisotropy 
and of x\\- For the small samples studied in this work, the size effects are still extremely 
important: for A^ = 9 the top is spherical, for A^ = 12 no definite sign can be ascribed to 
(27 — 17") • -^^^ ^^^ ^^'^ larger samples (A^ = 21, 27) there is a tendency towards a behavior 
of an oblate top (27 — 27~) ^ *-*' ^^^ ^^^ expected degeneracies of the symmetric top are not 
present leading to large uncertainties on x\\ (see fig-D . Unfortunately the Hilbert space of 
the quantum top in the A^ = 36 case is too large to allow the determination of the complete 
spectrum of low lying levels and we are thus unable to decide clearly even on the sign of the 
above quantity in this case. 

D. Spatial symmetries of the quantum effective top 

We have determined in subsection III.B the three IRs characterizing Neel order on the 
triangular lattice: 

El : [ k = 0, 7^^^ = ^, 7^2./3^ = ^, a,^ = ^ ], Fs : [ k = 0, 7^^^ = -^, 7^2V3^ = ^, 
a^^ = ^ ] and Fa : [ k = ±ko, 7^2V3^ = ^, (^x"^ = ^ ]• 

These IRs should be identified with the three IRs of the invariance group C^v of the mag- 
netic arrangement. Each symmetry of the lattice results in a permutation of the magnetic 
sublattices. Because the order parameter is an element of 50(3), such a transformation may 
be cleared by a global rotational symmetry of the spins. This property allows us to compute 
the number of occurrences of the three IRs (see Appendix B): 

11 



npi = (a + 36 + 2c)/6, n^^ = {a — 3b + 2c) /6, and rir^ = (a — c)/3, 

where 

a = 25 + 1, b = cos{S7i) and c = sin(— (25+ l))/sin( — ). (13) 

o o 

In the hypothesis of an a la Neel symmetry breaking, the low energy spectrum of Eq.|l] should 
thus contain for each S, Ms subspace: 

nn [ k = 0, 7^^* = ^, 7^2V3^ = ^, a,^ = ^ ] IR, 

nr, [ k = 0, 7^^^ = -^, 7^2V3^ = ^, o^x^ = ^ ] IR and 

npg couples of degenerate [ k = ±ko, 7^27r/3^ = "^ , cTx"^ = "^ ] IRs. 

This appears to be true in all the exact spectra that we have computed, for all S values 
up to A^/6 (see Table | and fig.^. 

The mapping of the low lying levels of Eq.^ on those of Eq.|T2] would imply a "quasi- 



degeneracy" of some Fi and r2 levels. This phenomenon is not present in the studied samples 
(see fig.^; this does not exclude the possibility of a quasi-degeneracy in the thermodynamic 
limit, but this explains the difficulty to extract x\\ from these data. 

At that point, we have determined both the dynamics and the symmetries of the QDJS 
that should appear in exact spectra of a system exhibiting a Neel ordered phase in the 
thermodynamic limit. As regards these criteria the spectra of the THA point in favor of a 
Neel ground-state: the scaling of Ij_ seems even to dismiss the case of quantum criticality 
(/_L should then scale as N^^'^ as shown by Azaria et al |4^). We shall come back to this 
point in the conclusion. 

E. Spectra of samples with N = 3p + I 

The periodic samples which do not frustrate Neel order, must have three sublattices 
invariant in the periodic boundary conditions, that is a number of spins A^ = 3p. On the 
nowadays computers and algorithms, memory requirements limit the studies to samples of 
A^ = 9, 12, 21, 27 and 36. In order to enlarge the number of data available and have a com- 
plementary check of our hypothesis we have thus relaxed the periodic boundary conditions 
to allow the studies of the 7, 13, 16, 19, 25 and 28 samples. Let us look at such a tiling of the 
infinite lattice (see fig.^j). In order to preserve the possibility of a three lattices symmetry 

12 



breaking on the infinite lattice, the translations operations Ti and T2 should be linked to a 
rotation of the spins: to be specific, in fig.^ the interaction between the spin (10) and (1') is 
chosen to be equal to the interaction of (10) with the periodic image of (1) rotated by 27r/3 
and so on. Such boundary conditions, which imply the selection of an axis of rotation, break 
the rotational symmetry of the hamiltonian and the total spin is no longer a good quantum 
number, but the component S3 of the total spin on the axis of rotation remains a conserved 
quantity. If we keep a global reference frame for the spins, the translational sjTiimetry is 
difficult to handle. But this essential symmetry becomes again trivial by the choice of local 
frames rotated by ±27r/3 (resp. =F27r/3) in one-step translations of ui (resp. U2). Within 
such a choice the Heisenberg hamiltonian reads 

n = 2Y^ s,.7^^^(27^/3)sJ7^3(27r/3), (14) 

where Sj, Sj are the spin operators in the local reference frames and 7^3(27r/3) = 
exp(— i7r/3cr3), era being the third Pauli matrix. We can in fact solve this problem for 
any value (0, ip) of the rotations associated to one step translations in the Ui and U2 di- 
rections (as long as the total rotation around a plaquette is zero). We have indeed done it 
and verified that the absolute minimum of the energy is obtained for the values (±27r/3) for 
samples with N = 3p + 1 and for and ±27r/3 for samples with N = 3p (see fig.P). The 
same method gives information on the spin stiffnesses of the system [^ . 

If the system has Neel order at T = 0, the spectrum of the low lying levels of Eq.[l^ 
is easily deduced from the general considerations of section III.C The present boundary 
conditions induce a cylindrical symmetry of the problem around the 3 axis of the laboratory: 
this implies that the helicity axis 3 of the Neel state (see fig.^ is now constrained to coincide 
with the 3 axis of the laboratory. There is only one degree of freedom left, associated with 
the rotations around this axis, and one constant of motion: the 6*3 component of the total 



spin. The effective hamiltonian Eq.ll in these conditions reduces to: 



02 c2 

the two other terms in Eq.[ll| not related to a constant of motion should average to zero 
in any eigenstates. In our problem, S^ can take all the values ranging from —N/2 up to 
N/2 and the degeneracy of each eigenstates of Eq.|l5| is 2 for 5*3 7^ and 1 for 5*3 = 0. 

13 



Diagonalization results entirely corroborate these deductions as can be seen in fig.|TO[ The 
analysis of the low lying levels obeying Eq.|T3| gives a determination of x\\ (see fig.|^). As for 
the N = 3p case, the values obtained by this method are in agreement and only a bit larger 
than the first-order spin-wave approximation, moreover the size effects seem to be roughly 
the same in the two approaches. 

F. Extrapolation of the ground-state energy per bond 

For ordered systems, the finite-size effect on the A^-ground-state is mainly due to the 
cut-off of the long wave-length excitations. In the supposed-to-be Neel order these exci- 
tations are magnons, linear in k and the leading term of the finite-size corrections to the 
ground-state energy is of order N~^^'^ and thus the ground-state energy per bond varies 
as N~^/'^ : more sophisticated approaches would allow to compute the first coefficients of 
these expansions as function of the spin- wave velocities p6|j48|j43|l . In fact in the previously 



developed theories [0,^] only the N = 3p samples have been taken into consideration: in 
this case the fluctuations associated to the three modes of magnons equally contribute to 
the renormalization. In the N = 3p + 1 samples, the approach to the singular points at 
k = and ±ko is different and we expect different coefficients in the A^~" expansion. It 
is effectively what is seen in the finite-size scaling analysis of the spin-wave susceptibilities 
(fig.|^) and of the spin- wave energies (fig.|ll|). 

The analysis of the diagonalization results is more subtle, because it depends on the 
parity of A^. The ground-state is either an S* = or an 5* = 1/2 state and in this latter case 
one has to take into account the top inertial effects described in the previous subsections. 
This leads to an S/N correction to the ground-state energy (an S/N^ correction to the energy 
per bond) which is noticeable for the small sizes considered here. In order to compare the 
exact results to finite-size spin-wave results, it is necessary to extrapolate the exact E{S, N) 
to an effective value E{0, N). Two estimates of E{0, N) can be computed. The first one is 
obtained by the subtraction from the ground-state energy of the main inertial contribution 
of < \l/o|7Yeff|^o > that is for N = 3p samples: 

< 2si.Sj >[Jlv=< ^o|2si.Sj|^o > - < ^ol^^^^^ — 1^0 > (16) 
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and for A^ = 3p + 1 samples 

C2 

< 2si.s,- >[Jlv=< ^o|2si.Sj|^o > - < ^ol^^^l^o > • (17) 

The second estimate is an average of the same quantity taken on all the QD JS up to S' = A^/6 

< 2s,;.s, >!,'1^= 1/Nqdjs Y. < QDJS\^^^\QDJS > (18) 

S=min 

where Ti is the Heisenberg hamiltonian Eq.|l] and Tiefr is given by Eq.l2 or Eq.l5 according 
to the number of spins in the samples. 

The first estimate (used in |3^) is different from the exact ground-state for odd N samples 
only, whereas the second estimate which is an average over a large number of levels, always 
differ from the exact ground-state. These averaged results are compared to finite-size spin- 



wave results (see fig.0) where it is seen that the exact results are lower than the spin-wave 
results, the finite-size effects on the two sets of data being nevertheless roughly the same. It 
has been shown by Azaria and coworkers that the finite-size correction to the ground-state 
energy of the N = 3p samples is [|^ : 

En = E^- a(c|| + 2cx)Ar-3/2 (19) 

where a is a geometrical form factor and cy , c± are in plane and out of plane spin- wave 
velocities. 

On the basis of our numerical data we cannot decide on the exact value of the renormal- 
ized spin-wave velocities: it does not seem to be largely different from the spin-wave results, 
both in the N = 3p (fig.[Tl|a) and in the A^ = 3p + 1 cases (fig.|lT]b). The complete set of 
data analyzed within the following hypothesis 

C\\ = C\\sw c± = C±sw (20) 

gives an estimate of the energy per bond in the thermodynamic limit 

< 2si.Sj >oo= -0.363 (21) 

(see fig.|ll]c). 

At the end of this section, we can conclude that the spectral properties of the Heisenberg 
hamiltonian give strong argument in favor of an ordered ground-state and that the numerical 
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spectral data and finite-size corrections are consistent with the picture and not very different 
from the spin-wave results. We will now discuss the important issue of the value of the order 
parameter of these Neel states. 

IV. ORDER PARAMETER OF THE NEEL STATES 

The check for long range order on the lattice could involve the measurement of the two- 
points correlations < Sj(0).Sj(r) > or that of a macroscopic observable as the sublattice 
magnetization. 

In view of the lattice sizes that can be studied, the analysis of the two-points correlation 
functions is rather unconclusive, the asymptotic behavior is far from reached and the sizes 
are still too small to check the Kennedy-Lieb-Shastry inequality [^. The measurement of 
the squared sublattice magnetization TVI^: 

M^ =< S\ >=< S| >=< S^ > (22) 

or of the static structure factor 5(k) defined in state |\E' > as 

5(k) =< ^|S_,..Sfc|^ > (23) 

(with Sfc in Eq.6) is easier. S(k) measured in the ground-state of the N = 3p samples 
is shown in fig.|l^. The signature of Neel order in the A^-ground-state should appear as a 
macroscopic value of M. or 5(ko) using the following identity: 

2 A^5(ko)+<^|S2|^>=9A^2 (24) 

It should be kept in mind that, in the studied samples, the largest sublattices spins 
range from 3/2 (for A^ = 9) to 6 (for A^ = 36): they cannot be treated as classical quantities. 
A4^ is the square of a vector, it should be normalized by its maximum quantum value 
(A^/6)(A^/6 + l) (reached in the fully aligned classical Neel state |^o(i, 5", Ms) > of subsection 
III. A). The structure factor (same tensorial form as Ai"^) should be identically normalized. 
Two definitions of the order parameter {X or Y) seem equally valid: 

36A^ y-=8^M, (25) 
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X saturates at one in the classical Neel state (and in the A^ = 9 sample where the quantum 
fluctuations are ineffective to reduce the magnetization sublattices) and should be decreased 
by the quantum fluctuations in the quantum ground-state []• For the N = 3p + l samples, it 
is impossible to define a sublattice magnetization. On the other hand one can measure the 
sum Ait of the magnetizations computed in the local basis: for a fully aligned Neel state 
the square of this vector should saturate to {N/2){N/2 + 1). We thus define in this case: 

X2 = f-^* ^ (26) 

N{N + 2) ^ ^ 

Finite-size scaling analysis indicates that the leading correction to these parameters 
should go as N~^^'^; analysis of the results (see fig.|l3]) shows that for these small values of 



A^ the subleading correction is important: the N ^ oo extrapolation is thus rather difficult, 
but an extrapolation to a zero value seems highly improbable. 

On the other hand, the comparison between the diagonalization results and first-order 
spin-wave results is interesting: it is seen in fig.|l^ that the two sets of results do not differ 
by large amounts. It should be noted that Miyake spin-wave results are in our point view 



incorrect for finite sizes H^. We agree with Deutscher and Everts formula for first-order 
spin-wave results which restrict the renormalization by quantum fiuctuations to k 7^ 0, ±ko 
wave vectors [^]. As it is well known the spin- wave hamiltonian cannot be bosonized for 
k = 0, ±ko (the Bogolioubov transformation becomes singular). It is a bit lengthy but 
straightforward to show that these three Fourier components of the spin-wave hamiltonian 



^Examination of this problem of normalization has been underestimated in previous studies 
1 27, 32 1 . Because of an erroneous normalization most of the weight of the extrapolation proce- 



dure has been carried on unphysical numbers larger than the saturation value of the fully aligned 
classical Neel state. This essentially explains the difference in the conclusions of the previous au- 
thors and ourselves as regard the question of order of the THA. Deutscher and Everts have been 
faced with the same kind of difficulty in their recent work on the Ji — J2 hamiltonian on the 
triangular lattice |^^; by studying the behavior of the system when J2/J1 -^ —00 they arrive at 
the same kind of renormalization to order 1/A^ than ours. The above line of reasoning is probably 
more general. 
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can be recombined to give the total spin S"^ and that they do not participate to the renor- 
mahzation of the order parameter (such a remark has recently been developed by Zhong and 
Sorella in the case of the square lattice Heisenberg hamiltonian ^l|]). Careful examination 
of these singular terms shows that the first-order correction to the ground-state energy is 
exactly obtained by the summation of the usual formula on the A^ points of the Brillouin 
zone: 

<2s,.Sj >,^=-l/4 + — ^(cJk-1) (27) 

k 

where cjk = (1 — 7k)^^^(l + 27k) ^^^ and 7k = /(k)/6 (see Eq.|^). On the other hand, 
the correction to the Aiz magnetization involves exclusively the k 7^ 0, ±ko points of the 
Brillouin zone that is A^ — 3 points in the N = 3p samples and A^ — 1 points in the N = 3p+l 
samples. 

1 1 1 -I- 2k 

<><«>= ill -]v 5: (^-1)1 (28) 

k^O.itko 

The careful comparison between the exact results and the first-order spin-wave results 
leads us to conclude that the spin-wave approximation seems an extremely good quantitative 
approximation for the considered sizes; on the basis of the present data, it seems highly 
hazardous to adopt other estimate of the thermodynamic limit that the spin-wave results 



^ □• In a previous paper |^3[,we have shown that, for a given size, the order parameter 



modulus is roughly the same in all the QDJS. 

V. CONCLUSION 

The thorough analysis of exact results of diagonalization of the Heisenberg hamiltonian 
on periodic samples with triangular symmetry has brought various pieces of information 



^Definitions and comparison of order parameter in symmetry breaking solutions (spin-wave ap- 
proximation) and symmetry non breaking one (exact diagonalizations) has given birth to many 
studies [p5tp4 -|55|. In this work, we take care to compare two parameters that are normalized to 
1 in the absence of fluctuations and should be equal in the thermodynamic limit. The problem of 
the symmetry breaking will be studied in a forthcoming paper. 



both qualitative and quantitative. 

The quahtative information emerge from the symmetry analysis of the spectra: this 
analysis shows a strict hierarchy of levels. The first family, degenerated to the absolute 
ground-state in the thermodynamic limit, contains all the quantum states and only those 
needed to build symmetry breaking Neel states with three sublattices magnetizations in a 
C^y invariant configuration. Their dynamics can roughly be mapped on the dynamics of 
a quantum symmetric top as expected from general considerations on the nature of the 
order parameter in the Neel ordered states Q. Understanding this family of levels from the 
coupling of three macroscopic spins explains why the multiplicity of this quasi-degenerate 
set of levels (the QDJS) is of order A^'^ (in the square lattice case the multiplicity should 
be of order N'^) which gives in the thermodynamic limit an entropy of order log A^. Rather 
unexpectedly the present analysis does not give quantitative answers extremely different 
from first-order spin-wave results. On the basis of the coherence of the data concerning the 
symmetries of the QDJS, their dynamics, the energies, the spin susceptibilities, and order 
parameter of both the A^ = 3p and N = "ip + 1 samples and their size effects, one can 
only sustain the hypothesis of an "ordered ground-state" for these small samples. Indeed we 
have no information of the effect of softer quantum fluctuations with wavelengths larger than 
about 6 lattice sizes. The clusters expansions of Singh and Huse and the high temperature 
expansions of Elstner, Singh and Young point to a weaker order (if any?) than obtained 
in the first-order spin-wave calculations. In these two last methods the invoked clusters are 
smaller than our largest sizes but there are no bias due to periodic boundary conditions and 
the sampling of the quantum fluctuations is different from what is done in this work. 

The present state of the art does not exclude that for larger sizes than those explored 
today quantum fluctuations drive the system towards criticality but there is in the small 
samples studied here no evidence of such a behavior. 



^We should mention that the first exhibition of such a structure of low lying levels invoked by 
Anderson in his seminal paper of 1952 appears in a paper by Gross, Sanchez- Velasco and Siggia 



1 34 1 that regrettedly escaped our vigilance in our preliminary report of 1992, indeed most of the 



theoretical hints were present in the first paper of this series. 
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VI. APPENDIX A: GROUP THEORY ON THE THA, NUMERICAL METHOD 

Here, we explain how to obtain the complete spectrum of the spin-1/2 THA, for the 
largest possible samples compatible with the triangular symmetry of the infinite lattice. 
The number of eigenstates increases as 2^ with N. Even for very small values of N, there is 
no stable numerical method to calculate directly all the eigenvalues of 7i. The problems to 
be handled are memory space requirement, computer time, and degeneracy of states. These 
problems can be reduced by an intensive use of space and spin symmetries. Indeed, when 
all symmetries are accounted for, one can work in subspaces where eigenvalues are no more 
degenerate and better control on round off errors is obtained. These subspaces come from 
the decomposition of the Hilbert space according to the irreducible representations (IRs) of 
the symmetry group of the problem. 

This symmetry group is the direct product of the space symmetry group of the lattice 
times the spin rotation group SU{2). Thus, its IRs are the tensor products of space-group 
IRs times SU{2) IRs. We first consider in details the space symmetry, then quickly the spin 
symmetries and explain how they are implemented in the numerical approach. 

The space symmetry group 

The spins stand on a triangular lattice defined by the basis vectors ui and U2 (see 
fig.§). On the infinite triangular lattice, the group Qoo of symmetries, which leave the 
lattice sites globally invariant and keep the vicinity relations, is the semi-direct product 
of the translation group 7^ times the point group V. The translation group T^ has two 
generators: the translations of vector Ui and U2. The point group V consists of the symmetry 
transformations keeping a site invariant: its generators are the planar rotations 7^^ and 7^27r/3 
and the axial refiection a^ with respect to ui (see fig.§|). The point group V is isomorphic 
to Cgt, (or "Dg) and has twelve elements: six rotations and six refiections. 
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In order to preserve the translation invariance and to reduce the number of degrees of 
freedom, periodic conditions are used, defined by the vectors Ti = /ui + mu2 and T2 = 
Z'ui + m'u2 (see fig-H). TZn is always a symmetry transformation of such a system. For 7?.27r/3 
to be a symmetry transformation of the periodic lattice, we chose T2 = T^tt/sTi (/' = — m; 
m' = I + m). The number N of sites per cell (Ti, T2) is thus given by: 

N = 1^ + lm + m^. (29) 

The system has an axial reflection when Im = ot I = m. Therefore the point group is 
V = Cq^ (resp. Cq if there is no a^ axis) and the space group Qn = TJv A Cq^ (resp. TJv A Cq) 
has 12N (resp. 6A^) elements. In the following, periodic samples with A^ multiple of 3 are 
used. 

Let E^ be the space of wavefunctions for a system of A^ spins. A wavefunction is a linear 
combination of configurations: 

^ = ^ac|c> (30) 

c 

where a configuration |c > is an element of the basis Bn = {|T>)|i>} of -^a^- -^a^ is a 
linear representation of Qj\f. It can be expressed as a direct sum of IRs of Qn- Summing the 
isomorphic IRs leads to a unique factorization in a direct sum of subspaces Ej^^p associated 
to the various IRs p of Qn- Each of these subspaces is invariant under Ti. Because Qn is not 
commutative, some of the IRs have dimensions greater than one resulting in a degeneracy 
of Ti on Efyf^p. Thanks to the rather simple structure of group Qn, one can eliminate this 
degeneracy by working on some subspace E'j^ associated to some one-dimensional IR p of 
some subgroup Qn,p of Qn- 

By using E'j^^ instead of E^^p, the eigenvalue multiplicities are divided by the dimension 
of p, which is also equal to the ratio of the cardinal of Qn to the cardinal of Qn,p- 

Let us first consider the case where there is no axial refiection: the point group V is an 
abelian cyclic group of six elements. Because 7/v is abelian, its IRs are one- dimensional and 
are characterized by k, a vector of the Brillouin zone. Thus, E^ is factorized in A^ subspaces 
ii^7v,k- The group V acts on these vectors to form different stars. Two vectors ki and k2 
belong to the same star if it exists an element of V (here a rotation) which transforms ki in 
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k2. Because this transformation commutes with the hamiltonian, the eigenvalue spectra of 
i?7v,ki and -EAr,k2 will be identical. For a fixed k, let V^ be the subgroup of V which keeps k 
invariant. This subgroup is abelian and all its IRs are one-dimensional. The cardinal of V^ 
is 6 for k = 0. It is 3 for k = ±ko (±ko are the corners of the Brillouin zone); this occurs 
only when N is multiple of 3. It is 2 if k = —k 7^ (the middle of the side of the Brillouin 
zone) and this occurs only when A^ is multiple of 4. It is 1 for the other k -vectors. The 
different IRs are listed in Table |I|. 

When there is an axial reflection {Im = 01 1 = m), the group V is generated by ax and 
TZn/s- If the group V^ is abelian, the previous construction is applied. If the group V^ is not 
abelian, some of its IRs are two-dimensional. This is the case for k = {V^ = V = Vq) and 
k = ko (Vk = V-s). The abelian normal subgroups are respectively Cq and C3. The different 



IRs are listed in Table III 



For small systems, other accidental space symmetries can occur which imply new de- 
generacies. If they commute with all the previous space symmetries, the degenerate states 
appear in the same IR. Otherwhise, the degenerate states stand in different IRs. 

The spin symmetry group 

The hamiltonian of Eq.|l] is invariant under global spin rotations: it commutes with S^, 
where S is the total spin. For spin-1/2, the symmetry group is SU{2). The IRs are labeled 

hjS. 

Numerical method 

A general wavefunction is defined by the knowledge of the set of 2^ coefficients {ad 
Eq.|30|. The first easy reduction of this set is to work in an eigenspace of S^ '■ the basis 
size thus becomes C^ " . The next step is to work in a given E'^ , characterized by the 
phase factors associated to each element of Qn,p- One groups together the configurations 
in conjugate classes: \c > and \c' > belong to the same class if it exists an element X of 
Qn,p such that \c' >= X\c >. So, Oc' differs from ac by a phase factor listed in Tables |I| 
or |T|. The number of independent coefficients is therefore reduced by roughly a factor A^ 
(translations) and at maximum by a factor 12A^ (cardinal of Qn)- The memory requirement 
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and computer time are reduced by the same amount. 

Implementation of the SU{2) symmetry is done via a projector technique by using the 
operator: 

A A S{S + I) -^i^ + l) ^^^^ 

i=N/2,N/2-l,... ^ / \ / 

where S^ is computed thanks to relation S^ = S'^ + Sz + 2S-S+ and S+ and S- are computed 
like 7i. Applying S^ is therefore as fast as applying the hamiltonian, or so. 

Lanczos method has been applied to diagonalize the hamiltonian in the subspace asso- 
ciated to each IR p. Working with classes instead of configurations insures that the vectors 
generated by Lanczos method stay in the same subspace. For small systems, the hamiltonian 
and the two operators S^ and S^ are tabulated; for large systems (A^ = 36), they have to be 
computed at each step, but vector components are stored in an order allowing application 
of these three operators to be vectorized. In any case, diagonal part of the hamiltonian is 
computed separately at once and stored. Additional details on the technical tricks will be 
given elsewhere. 

The algorithm is as follows: i) first, for some IR, build the conjugate-class table and 
the phase-factor table; ii) tabulate 7i, S+ and S-; it) choose a random initial vector and 
project it out in the desired subspace of S^; iv) apply the Lanczos algorithm; v) calculate 
the eigenvectors, correlations... 

In fact, round-off errors propagate very quickly in this algorithm. In order to eliminate 
these errors, computed vectors are orthogonalized to all previous ones and projected out in 
the S'-subspace. One can compute the dimension rriplS) of this subspace and verify that, 
after exactly mp{S) iterations, the modulus of the last vector is zero. This is a strong test 
of this algorithm. 

VII. APPENDIX B 

Here, we compute the number of replica of each irreducible representation of the QDJ 
states if a Neel state occurs. The idea is to use some symmetries that keep invariant the 
classical Neel state. These symmetries are the compositions of the permutations of the 
sublattices times specific spin rotations. 
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Thus, let us start with an a priori classical Neel state and its symmetry group. In our 
case, the Neel state is made of three sublattices as described in section II. Let us define the 
plan (xy) of the classical Neel state and z a perpendicular axis. The group which permutes 
the three sublattices is 5*3 (isomorphic to the dihedral group D^). After a permutation 
which exchanges the spins between the sublattices (giving an other Neel state), the initial 
Neel state is recovered if the spins are rotated simultaneously. For example, the permutation 
of two sublattices (say A and B) times a spin rotation of axis parallel to the spins of the 
last sublattice (thus C) recovers the initial Neel state. Likewise, a cyclic permutation of the 
three sublattices times a spin rotation of 27r/3 along the z-axis also recovers the initial Neel 
state. Thus, two (resp. three) sublattices permutations are associated to an half (resp. a 
third of a) turn around an axis in the (xy) plane (resp. around the 2;-axis). 

In the quantum case, the permutation of the sublattice is one of the symmetry of the 
lattice (say, central symmetry for the two sublattice permutation and translation of one step 
for three sublattice permutation) times a spin rotation in SU{2). The group of these spin 
rotations is the dihedral group D3 for integer spins {SO (3)), and is a group of 12 elements 



for half-integer spins {SU{2)), whose character table is given in Table |rv|. Note that the type 
of spin rotation associated to each class is specific of the Neel state and is necessary (here 
also sufficient) to calculate the number of replicas rir^ of each irreducible representation (IR) 
Fj. As all the above transformations conserve the total spin S, one can thus calculate the 
traces of these rotations in the subspace of fixed S, Ms: 



^"' = ^ 5Z^r(^)(7^fc)x.(A;)A^,,(A;) (32) 



k 



For integer (resp. half-integer) spins, only the nr^ (resp. rap') are non zero. 

Now, we identify the Fj IRs with the IRs of the hamiltonian. Because the QDJS have 
already been found to belong to three precise IRs (see section III.B), this identification is 
straightforward: Fq is the trivial IR, Fi is the IR odd with respect to the inversion (k = 0, 
invariant under a 27r/3 rotation), and F3 is the last doubly degenerate representation (k = ko, 
invariant under a 27r/3 rotation). The identification of the F^ is the same for half- integer 
spins and leads to formulas (13) both for integer and half-integer spins. 
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TABLE I. Lowest energies, degeneracy and quantum numbers for the samples A^ = 9 , 12, 21, 
27, 36. Components of vectors k are in units of 2tt/N. In the 3 last columns, 1 means invariant 
under the symmetry, means no symmetry and -1 means a phase factor under the symmetry (j for 
the rotation of 27r/3 and -1 for the two others). Stars stand for the Quasi-Degenerate- Joint-States 
(QDJS). 
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^27r/3^=J^ 
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< ^27r/3 > 


^27r/3V'=V' 
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^27r/3^=J^ 
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k = -k / 


{id,n^} 


7^^V'=±V' 


3 


other k 


{Id} 


6 



TABLE IL Irreducible Representations (IR) when there is no axial reflection {Im j^ or I j^ m). 
First column: vector of the Brillouin zone. Second column: subgroup which keeps this vector 
invariant; < 7^27r/3 > stands for the group generated by TZ2-n/z- Third column: IR list of this 
subgroup, with the phase factor associated to the transformation; ^ is a wavefunction. Fourth 
column: multiplicity of each eigenvalue found in this IR. 
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6 
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TABLE III. Idem as in Table ||, when there is an axial reflection {Im = or / 



m . 
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TT 
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d 
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sm(fd) 
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sm(f(i) 



TABLE IV. Character table for the spin-rotation group of SU{2) involved in the Neel state; p2 
and p3 denote the class of the spins rotations associated to 2 and 3 sublattice permutations. The 
number of elements in each class is N^i. The rotation angle associated to each class is </> and the 
trace of the rotation in the subspace of fixed spin S is given in the last line, where d = 25* + 1. 
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FIGURES 
FIG. 1. The classical Neel ground-state: on each ijk triangle Si + Sj + s^ = 0. This defines 
three sublattices A, B, C on which the spins are ferromagnetically aligned; the angle between the 
spins of two sublattices is 27r/3. For a given planar upwards triangular plaquette described in the 
counter-clockwise direction, the spins can rotate clockwise or counter-clockwise, corresponding to 
the two different helicities: here a positive helicity is assumed. 

FIG. 2. Energy spectra of Eq.|| versus S^ = S{S + 1) ; a) complete spectrum for A'^ = 9 and 12; 
b) lowest energies for A^ = 21 and N = 27. The horizontal and vertical scales have been enlarged 
by the same factor so that the slope of the energy per bond versus S(S-|-1) can be compared. One 
sees on these graphs that this slope goes rapidly to zero. The straight line is a guide for the eye to 
link the low lying energy levels called QDJS (for quasi-degenerate joint states). 

FIG. 3. Moment of inertia versus sample size 

FIG. 4. a) Enlargement of the low energy levels for A^ = 9 versus S{S + 1). Three-legs-symbols 
family and - family are essential ingredients of an a la Neel symmetry breaking. They respectively 
span the subspace of quasi-classical Neel ground-states and the subspace of the long wave-length 
k 7^ excitations: i.e. the magnons. Black triangles represent the chiral states. Open triangles 
are states invariant under rotation of 27r/3 and odd under axial symmetry, b) same as in a) for 27 
spins, c) same as b) with ^^^j ' subtracted from the whole spectrum. 

FIG. 5. Schematic of the free dynamics of a symmetric top 

FIG. 6. Comparison between the exact spectrum and a fitted symmetric top for the N = 27 
sample : 1, 2 and 3 stand for the Fi, r2, T^ IRs (see text); □ : spectrum of an ideal symmetric top 
(the isotropic term }j^j ' is subtracted from the spectrum to focus the comparison on the second 
term of Eq.l2); For N = 27, all levels of the ideal top (□) are doubly degenerate (we do not take 
into account the trivial 25* -|- 1 magnetic degeneracy). In the exact spectrum, the levels 1 and 2 
are simply degenerate and the 3 levels are doubly degenerate. Perfect agreement between the two 
spectra would necessitate a quasi degeneracy of 1 and 2 levels. 
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FIG. 7. Spin susceptibilities of the THA on finite samples. Fig. 7a: xi. (triangles) and x\\ 
(square) normalized by its classical value for the N = 3p-samples; large symbols stand for results 
obtained from the QDJS analysis; small symbols show the finite-size results of spin-wave calcu- 
lation of Chubukov et al [^]. Points represent spin- wave results and indicate the infinite size 
extrapolation. Fig. 7b: x\\ of the N = 2>p + 1-samples; symbols are the same as in a). 

FIG. 8. The A^ = 13 sample and the tiling of the infinite lattice. Numbers stand for the spin. 
Prime (resp. double prime) means that the original spin is rotated by ib27r/3 (resp. ib47r/3). 

FIG. 9. Opposite of the ground-state energy of the A'^ = 21 and A^ = 19 samples as a function 
of the spin-rotation angles 0, V' attached to the translations ui, U2. The spectrum depends only 
on the angles ^ = 1(f) + mip, ^ = —mcf) -|- (/ + rn)%l) attached to the translations Ti, T2. In the 
N = 3p case, the absolute minimum of the energy is obtained for angles <I> ,^ equal to or :iz-^ 
characterizing the classical Neel states and only (<5, ^) = (±-|-, -F"f") in the A^ = 3p + 1 samples. 

FIG. 10. Spectra of A^ = 3p -|- 1 samples versus S'^. Note the tower of states, only doubly 
degenerate, and its collapse to the absolute ground-state with increasing A^. 

FIG. 11. Finite-size scaling of the ground-state energy per bond and comparison with spin-wave 
results, a) three-legs symbols: diagonalization results for the A^ = 3p-samples; black triangles: 
averaged values computed thanks to Eq.l8; small open triangles: spin-wave results; dotted line: 
7\r-3/2 £^ Qj-^ large A^ spin- wave results, b) same as in a) except crosses: diagonalization results for 
the N = 3p + 1-samples. c) averaged values of < 2sj.Sj > versus N~^''^; dotted line : N^'^'"^ fit to 
these results of diagonalizations (see text). 

FIG. 12. Static structure function in the ground-state of the N = 3p samples versus the k- vector 
modulus. 

FIG. 13. Neel order parameter as a function of the sample size. The order parameter is nor- 



malized by its maximum value (see Eq.25 and Eq.2£). Triangles (resp. squares) stand for A^ = 3p- 



(resp. N = 3p + l)-samples; black symbols stand for diagonalization results; open symbols stand 



for first order spin- wave results (2 < M^ > of Eq.28); dotted lines : large A^ fits (A^ > 5000) of 
first-order spin- wave results {X^ = X^o + aN^^'"^). 
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